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A parameter whose coupling to a quantum probe of n constituents includes all two-body inter- 
actions between the constituents can be measured with an uncertainty that scales as Xjr?^"^ , even 
when the constituents are initially unentangled. We devise a protocol that achieves the scal- 
ing without generating any entanglement among the constituents, and we suggest that the protocol 
might be implemented in a two-component Bose-Einstein condensate. 
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Quantum mechanics determines the fundamental lim- 
its on measurement precision. In the prototypal quantum 
metrology scheme, the value of a parameter is imprinted 
on a quantum probe through an interaction in which the 
parameter appears as a coupling constant [l[ . The num- 
ber n of constituents in the probe is often considered to 
be the most important resource for such schemes. We 
denote the parameter to be estimated by 7, and we write 
the interaction Hamiltonian as Ti. — h'jH , where H is a 
dimensionless coupling Hamiltonian. The measurement 
precision is quantified by the units-corrected root-mean- 
square deviation of the estimate 7est from its true value. 

The essential point made in 0, S i, H, H, 0, 11 H is that 
the scaling of ^7 with n depends on the probe dynamics 
as expressed in H. For an interaction that acts inde- 
pendently on the probe constituents, the optimal mea- 
surement precision scales as 1/n, a scaling often called 
the "Heisenberg limit," as this was believed to be the 
best scaling allowed by the Heisenberg uncertainty prin- 
ciple. In contrast, a nonlinear Hamiltonian that includes 
all possible fc-body couplings gives an optimal sensitivity 
that scales as l/n'"'. To achieve this requires that the ini- 
tial probe state be entangled. If practical considerations 
preclude initializing the probe in an entangled state, sen- 
sitivity that scales as l/n*^~^/^ is possible using a probe 
that is initially in a product state [3, H, 0, 0, @1 ■ Both of 
these scalings can be achieved with separable measure- 
ments. 

Practical interest in using nonlinear interactions for 
quantum metrology comes from the fact that even with 
two-body couplings and initial product states, it ispos- 
sible to obtain a l/n^/^ scaling for ^7 [2,, A M, 0, S H- 
In some versions of this scheme, entanglement is gener- 
ated during the protocol that leads to the l/rv^^^ scaling. 
We formulate here a protocol that generates no entangle- 
ment among the probe constituents, yet still achieves the 
Xjr?^"^ scaling; in this protocol, it is clearly the dynamics 



alone that leads to improvement over the 1/n scaling. 
Even though this Letter is mainly about improving on 
the Heisenberg scaling, any experimental demonstration 
of a scaling better than l/n-^^^ would be of considerable 
interest to the metrology community. 

A typical k = 2 choice for a probe made of qubits is 
H — , where Jz — \ ^3 is the z component of the 

"total angular momentum," with Zj being the Pauli Z 
operator for the j th qubit. If we denote the eigenvectors 
of Z by |0) and an optimal initial state is the en- 
tangled "half-cat" state, (|0 . . . 0) -I- 11 . . . 1))|0 . . . 0))/\/2, 
where the first part of the state refers to the first half of 
the qubits and the second part to the rest. This state 
evolves to (e-'T*"'/''|0 + 1))|0 . . . 0))/n/2 after 
time t. Measuring the product of Pauli X operators on 
the first half of the qubits gives a signal that oscillates in 
7 with frequency An estimate of 7 based on sam- 

pling from this signal over v trials leads to a measurement 
precision ^7 = 'i/tn'^^ 

If the initial probe state is required to be a product 
state, an optimal input state is of the form e^^'^a'^jO)®" = 
[cos(/3/2)|0) + sin(/3/2)|l)]®", where < ^ < 7r/2. A 
measurement of Jy after the probe has evolved for a time 
t under the Hamiltonian leads to a measurement pre- 
cision that scales as Xjtrc'^'^^/v, provided is small J^. 
This scaling applies for all values of /9 ^ 7r/2, but the 
optimal sensitivity occurs for fi = 7r/4. The restriction 
to small times arises because the Jz eigenstates in an 
expansion of the evolving state accumulate phase shifts 
quadratic in n, leading to a "phase dispersion" that after 
a short time renders it impossible to determine 7 opti- 
mally from a separable measurement such as that of Jy. 

The Jj? Hamiltonian is entangling. The entangle- 
ment generated during evolution from an initial prod- 
uct state and the phase dispersion are two aspects of the 
same phenomenon. One might think that the generated 
entanglement and associated phase dispersion somehow 
play a role in the enhanced scaling, but it would 

normally be expected that the phase dispersion is best 
avoided 

The essential observation we make here is that if the 
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Hamiltonian were replaced with one of the form H = nJz , 
there would be no phase dispersion and no generated en- 
tanglement. An nJz Hamiltonian acts as a linear cou- 
pling whose strength is proportional to n. Physically, an 
nJz coupling cannot arise from a fundamentally linear 
coupling, as that would require the coupling strength to 
be a function of the number of constituents in the probe, 
but it can arise naturally from quadratic couplings to the 
parameter. 

With a pure nJz interaction, the optimal initial prod- 
uct state is e~''^y''/'^\0)'»"- = [(|0) + |l))/y2]®". The 
state remains unentangled at all times, evolving to 
[(-g-«7t"/2|o) + e''^'*"/^\l))/V2]'^". A measurement of 
at time t has expectation value (Jx) — cos ^tn and un- 
certainty AJj; = sin7in|, leading to a measurement 
precision 6^ — /^.J^/ \/v\d{Jx) I d-y] = l/tn^/'^y/v after v 
trials. A measurement of any other equatorial component 
of J achieves the same sensitivity. The enhanced scaling 
in a protocol that uses an nJz coupling and an initial 
product state is clearly due to the dynamics alone, not 
to entanglement of the constituent qubits. These results 
indicate that in quantum metrology, entanglement is im- 
portant only in providing an optimal initial state, which 
leads to an improvement by a factor of over initial 
product states. 

We are interested in investigating measurement proto- 
cols that use both and nJz interactions in systems 
of bosons that can occupy two modes with creation op- 
erators a\ and a\. In the Schwinger representation, we 
have Jz — \{ni — 712) and n = ni-\- 712, where ni = a\ai 
and 712 = 12 "^2 ^'^^ numbers of particles in the two 
modes. The bosons we consider interact with one an- 
other, but the interactions conserve particle number, so 
the system has a nonzero chemical potential. Our mea- 
surement protocols, for both types of coupling, can be 
represented in terms of the interferometer with nonlinear 
phase shifters depicted in Fig. [T] In practical implemen- 
tations, the interferometer might be an optical or Ramsey 
interferometer or an interferometer made up of coupled 
nanomechanical resonators Q . 

An nJz coupling acts as a linear coupling with a cou- 
pling strength proportional to n. Thus the effect of de- 
coherence on our measurement protocol is the same as 
that on a linear protocol with a product-state input. In 
particular, decoherence that acts independently on the 
probe particles does not change the 1 /n'^/^ scaling 0, Q . 

We turn now to the problem of implementing the non- 
linear interferometer of Fig. [1] in a laboratory system of 
considerable interest. For this purpose [l, [sl we consider 
a two-mode Bose-Einstein condensate (BEG) in which 
the n atoms can occupy two internal states (modes) la- 
beled |1) and |2), which are typically hyperfine levels. 
The atoms that form the initial BEG are all in the in- 
ternal state In the mean- field approximation, they 
all share the same spatial wave function 'ipn[r), which 
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FIG. 1: Nonlinear interferometer giving jf and nj^ cou- 
plings. An incoming beam of n bosons is split at a beamsplit- 
ter, which puts each boson into an appropriate superposition 
of being in the two arms (modes). The two initial nonlin- 
ear phase shifters produce Kerr phase shifts Xi*^! s-iid Xaf^-i- 
The phase shifter at the intersection of the beams produces 
a cross-Kerr phase shift 2xi2'ii7i2. The final 50/50 beam- 
splitter converts the required measurement of an equatorial 
component of J into a measurement of Jz, i.e., a counting 
of the difference of the numbers of particles in the two out- 
put beams. The net effect of the nonlinear phase shifters is 
the same as a probe Hamiltonian Ti acting for a time t, with 
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Htlh = xiJi? + Xin\ + 2xi2nin2 

= (X + Xi2)nV2 + (Xi - X2)n3z + 2(x - X12) 

where x = ^(Xi + X2) is the average Kerr phase shift. The 
first term in the second line produces an overall phase shift 
and can be ignored. The nJz coupling comes from having 
different Kerr phase shifters in the two arms; to eliminate the 
interaction requires a cross-Kerr coupling xi2 = X- Under 
these circumstances, we have Ti = hrfnJz , with 7t = Xi ~ X2 • 
The case X2 = — Xi yields a pure nJz couphng without a 
compensating cross-Kerr phase shift. 



is the n-dependent ground-state solution of the Gross- 
Pitaevskii equation for a trapping potential V{r) and a 
scattering term characteristic of internal state An 
external field, playing the role of the first beamsplitter 
in Fig.^n drives transitions between the two internal 
states [10|, resulting in every atom being in the same 
superposition of the two internal states. We assume 
that the atomic collisions are elastic, so the only scat- 
tering channels are |1)|1) ^ |1)|1), |2)|2) ^ |2)|2), and 
|1)|2) |1)|2). These have amplitudes gn, 922, and 512, 
where gij — Anh^aij /m, with being the s-wave scat- 
tering length. The effect we seek is the differential phase 
shift between the two internal states due to their different 
scattering properties. After some period of evolution, a 
second external field, playing the role of the second beam- 
splitter in Fig.[Tl drives a 7r/2 pulse between the internal 
states. A final measurement then determines the popu- 
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lation difference between the two internal states. In the 
foUowing we are interested in the BEC dynamics that 
occurs between apphcation of the external fields. 

We assume that the two internal states are chosen 
so that both see the same trapping potential V, which 
is a situation that can be achieved in the laboratory. 
Nonetheless, the spatial wave functions corresponding to 
the two internal states will diverge because they experi- 
ence different scattering interactions. The effect of the 
scattering terms on the spatial wave functions becomes 
important at the atom number Uc where the scattering 
energy becomes comparable to the total atomic kinetic 
energy. For n small compared to ric, the two spatial wave 
functions remain essentially the same, and for n much 
larger than ric, the spatial changes, though they become 
substantial, occur on a time scale longer than the phase 
shifts of interest by a fractional power of n/ric, which can 
be around ten in laboratory experiments We thus 

neglect changes in the spatial wave functions, assuming 
that both internal states retain the initial wave function 
ipnir) for the duration of our proposed experiment. 

With these assumptions the Hamiltonian for the two- 
mode BEC [11, [11 takes the form 
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where ?] ^ J dr |^„(r)|'*, 71 = ^(gii -322), 72 =9-912, 
and g — ■^{gn +922) (notice that 71 and 72 do not have 
units of frequency). The only effect of the Hamiltonian 
Ho = nEo + j{g + gi2)rin'^ — ■^grjn, where Eq is the single- 
particle kinetic plus trap potential energy corresponding 
to ipn, is to introduce an overall phase, and thus 7io can 
be ignored. We assume n is large enough that we can 
replace n — 1 with n in Ti.. 

In a harmonically trapped BEC, the repulsive scat- 
tering interactions cause the single-particle ground-state 
wave function i/'n to spread as the number of particles in- 
creases. This effect appears in the BEC Hamiltonian in 
the factor 77, which is inversely proportional to the effec- 
tive volume occupied by the ground-state wave function. 
The n dependence of 77 gives the coupling strength a de- 
pendence on n that must be included in our analysis of 
the precision in estimating 71,2- 

When the number of atoms is small compared to ?ic, 
the total kinetic energy far exceeds the scattering en- 
ergy, resulting in a ground-state wave function that is 
independent of n. In a three-dimensional harmonic trap 
with ground-state half-width s, the total kinetic energy is 
^ n{h'^/ms'^), and the scattering energy is ^ n'^{gii/s^) 
(for atoms in internal state giving Uc ~ s/an. Typ- 
ical values of an ~ 10 nm and s ~ 10 ^ui give nc ~ 1 000. 
Hence, for a condensate composed of tens to a few hun- 
dred or so atoms, ij does not depend significantly on n, 
implying a scaling of Xjr?^'^ in such small BECs. 

In large harmonically trapped BECs, with n 3> Tic, '^ 
acquires an n dependence that defeats the desire to im- 
prove on 1/n scaling. Strategies for dealing with this 



include using traps with harder walls than a harmonic 
trap and working with BECs confined to fewer than three 
dimensions. To assess these strategies, we compute the 
n dependence of r] when the BEC is trapped in d lon- 
gitudinal dimensions by a spherically symmetric poten- 
tial V = and is tightly confined in the remaining 
D — 3 — d transverse dimensions by a harmonic potential. 
The longitudinal trap is characterized by the hardness 
parameter q and the half-width of its (bare) ground-state 
wave function, Rq = (ft^/mfc)^/^'^"'"^-', for which a typical 
value might be Rq 10 /xm. The tight transverse po- 
tential is characterized by its resonant frequency ujq and 
the half- width s — (fi,/2mijJo)^/^ of its ground-state wave 
function, for which a typical value for a tight trap would 
be s ~ 100 mn. 

There are now two critical atom numbers. The first, 
ul = {Rq / aii){s / Rq)^ , occurs when the scattering en- 
ergy is comparable to the longitudinal kinetic energy. 
As n increases from n^, the ground-state wave func- 
tion spreads in the longitudinal dimensions, its size grow- 
ing as i? ^ Ro{n/nLY/^'^^'^\ The second critical atom 
number, ut = (s/aii)(i?o/s)'^*-'~''^''^'', arises when the 
scattering energy becomes as large as the transverse ki- 
netic energy (and thus does not apply when d = 3), at 
which point the longitudinal extent of the wave func- 
tion is Rt- The corresponding atomic number density, 
utIs^R^ ^ 1/aiis^ ~ 10^^ cm"-^, is somewhat above 
the upper limit on number density set by three-body scat- 
tering losses. Thus we need only consider atom numbers 
smaller than n^. 

For atom numbers between and nr, a reasonable 
approximation to the ground-state solution of the Gross- 
Pitaevskii equation is obtained by using the Gaussian 
ground state of width s in the transverse dimensions and 
using the Thomas-Fermi approximation for the longitu- 
dinal wave function [3]. In this approximation we find 
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where X — jj, — ^Dhuio is the longitudinal part of the 
chemical potential /i and a^.d is a geometric factor of 
order unity that depends on d and q, but not on n. The 
n dependence of 77 implies an effective coupling strength 
that scales as 1/2, where C = (d + 3g)/2(d + g). The 
precision of estimating 71 or 72 thus scales as 

For a three-dimensional BEC trapped in a harmonic 
potential, the measurement precision scales as l/n^^^", 
worse than the Heisenberg scaling, but still better than 
l/ri^/2. To achieve super-Heisenberg scalings requires a 
trapping potential that is harder than a harmonic po- 
tential or else working with a one- or two-dimensional 
BEC. For d — 2, a. BEC trapped in a harmonic po- 
tential matches the 1/n scaling, and a one-dimensional 
harmonic BEC betters it, achieving a l/n7^^ scaling. 
A d-dimensional BEC achieves super-Heisenberg scaling 
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when the hardness parameter q exceeds d. The hmit of 
large q corresponds to a trap with hard walls and ex- 
tent 2i?o and has ^ = 3/2 regardless of d. For a one- 
dimensional BEC, an alternative to hard caps is to use a 
ring geometry. 

A good candidate for implementing the generalized 
metrology protocol is a BEC made of ^^Rb atoms. Atoms 
in the hyperfine level \F = — —1) = |1) are 

trapped and cooled to form a BEC, and then a Ra- 
man or microwave-driven transition is used to create 
a superposition of |1) and the hyperfine level [3, 15 1 
\F = 2;Mf = 1) = |2). The s-wave scattering lengths 
ail, 022, and ai2 are nearly degenerate for *TR,b, with 
ratios {022 : ai2 : an} = {0.97 : 1 : 1.03}. These val- 
ues imply that 72 = 5(511 + 522) — 512 is essentially zero 
for this scheme, meaning that a *^Rb BEC can realize 
the generalized quantum metrology protocol with a pure 
coupling. The optimal initial state for this pro- 
tocol has all atoms in an equally weighted superposition 
of |1) and 1 2). The quantity that is estimated is pro- 
portional to 71 = 5(511 — 522), which, though small, is 
nonzero for the scattering lengths in ^^Rb. 

Loss of atoms from the trap is an important decoher- 
ence mechanism, mainly due in our protocol to inelas- 
tic spin-exchange collisions (exchange of atoms with the 
thermal cloud that is present around any realistic BEC 
is negligible and can be ignored). A chief advantage of 
using protocols that do not rely on entanglement is that 
loss of atoms does not affect the sensitivity scaling, al- 
though it does generally degrade the sensitivity. In the 
case of spin-exchange collisions, the decoherence can be 
modeled in terms of a parameter Tri/2, which we can 
estimate using data from ll| and the assumption that 
|1) and |2) have the same spatial wave function. This 
estimate gives r/271 ~ 1/26, implying that we can per- 
form a measurement of 71 before inelastic collisions have 
a significant impact. 

A final issue is that the number of atoms in a BEC 
is not known to arbitrary precision, as we have assumed 
up till now. We propose to determine n by counting the 
number of atoms in both internal states at the output 
of our protocol. A determination of n with a fractional 
error of I^njn ^ 0.01, which is within current capabili- 
ties, would be sufficient for the purpose of demonstrating 
an enhanced scaling with n, provided the measurement 
time is kept short enough that the nonlinear phase shift 
is much smaller than n/An. We note that if An is big- 
ger than y^, the chief practical advantage of the nJz 
interaction is obviated, since the requirement on mea- 
surement time is as strict as or stricter than that set by 
phase dispersion in a protocol. Even so, the ability 
of the n Jz coupling to achieve enhanced scalings with no 
generated entanglement remains an important theoret- 
ical objective. Moreover, decoherence is likely to limit 
the measurement time more severely than either phase 



dispersion or number uncertainty. 

We have shown that it is possible to achieve mea- 
surement precision that scales better than I/71 by us- 
ing the dynamics generated by nonlinear Hamiltonians. 
The pure nJ^ scheme introduced here does not use quan- 
tum entanglement at any stage to achieve the enhanced 
scaling. Early experiments to test our scheme in BECs 
are likely to focus on demonstrating enhanced scaling in 
the estimation of some combination of atomic scatter- 
ing lengths. To be useful, however, our scheme must be 
adapted to measuring external fields that modulate the 
atomic scattering properties. One possibility is to use a 
i33Cs BEC with optical trapping of the = 3; Mf = 3) 
state [16J , which has a very broad Feshbach resonance at 
8 G, which makes the scattering lengths very sensitive to 

This 
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the strength of an external magnetic field 
suggests that our scheme might be used for ultra-high 
precision magnetometry. 
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